Temperature dependence of rho- and ai-meson masses 
and mixing of vector and axial-vector correlators 
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Within a chiral model which provides a good description of the properties of p and ai mesons in 
vacuum, it is shown that, to order T 2 , the p- and ai-meson masses remain constant in the chiral 
limit, even if at tree level they are proportional to the chiral condensate, ero. Numerically, the 
temperature dependence of the masses turns out to be small also for realistic parameter sets and 
high temperatures. The weak temperature dependence of the masses is consistent with the Eletsky- 
Ioffe mixing theorem, and traces of mixing effects can be seen in the spectral function of the vector 
correlator at finite temperature. 
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It is commonly believed that the enhancement of the 
dilepton spectrum in the invariant mass range around 400 
MeV observed in heavy-ion experiments Q is a conse- 
quence of medium modifications of the p-meson spectral 
function, i.e., more precisely, the spectral function of the 
vector-isovector current-current correlation function, in 
hot and dense matter. Most of the more "conventional" 
calculations find a strongly increasing p meson width (see 
Ref. H and references therein), but the more speculative 
hypothesis of a "dropping p mass" by Brown and Rho || 
is being discussed as a viable alternative 0] . The aim of 
the present study is to combine the phenomenology in- 
cluded in the more "conventional" calculations with the 
constraints imposed by chiral symmetry. As a first step, 
we have constructed in Ref. H a chiral model which gives 
a good description of the vector and axial- vector spectral 
functions in vacuum. The model is based on the linear 
a model, which has been extended to incorporate the p 
meson and its chiral partner, the ai(1260), as elemen- 
tary fields. In this letter we will discuss results from the 
application of this model to finite temperatures. 

Within the model presented in Ref. p] the vector me- 
son masses at tree level are given by 



+ h 2 cr 2 , m 2 = ml + (hi + h 2 )al 



(1) 



where mo is the vector-meson mass parameter in the La- 
grangian, hi and h 2 are coupling constants, and <jq is the 
expectation value of the scalar field a at tree level. As 
shown in Ref. ]5), the parameter too can take any value 
below w 400 MeV, in particular, to = 0. Thus, if Eq. (Q) 
were valid in general, i.e., not only at tree level, "Brown- 
Rho scaling" would emerge, i.e., the vector meson masses 
would go to zero at the chiral phase transition, where (a) 
vanishes. 

A good description of the experimentally measured 
vector and axial-vector spectral functions in vacuum 
can be achieved by including the one-loop p and a\ self- 
energies, which comprise the dominant decay channels 
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FIG. 1: Diagrams generating the dominant medium contri- 
butions to (a) (a), E p (b-d) and E ai (e-g). 



p — > 7T7T and a\ — > np as well as many other one-loop 
graphs which are necessary to comply with chiral sym- 
metry j^] . The divergences generated by the loops can be 
canceled by chirally symmetric counterterms, which are 
adjusted such that masses at one loop are equal to those 
at tree level ||. However, when we calculate the same 
diagrams at finite temperature T (i.e., when we replace 
the loop integrals over by a summation over discrete 
Matsubara frequencies 2-nnT\ the masses at one loop 
become temperature dependent, while the masses at tree 
level remain constant. 

Let us start by inspecting the temperature dependence 
of the condensate (a) . In the one- loop approximation it 
can be written as 



(a) = (T + Acr . 



(2) 



where Atro denotes the sum of all one-loop tadpole 
graphs, divided by im 2 ||. It can be decomposed into 
vacuum and medium contributions, Atro = Atro vac + 
AtTomed- The dominant medium contribution is gener- 
ated by diagram (a) shown in Fig. |l|. In the chiral limit, 
TOjr — > 0, it is of the order T 2 , while all other one-loop 
diagrams are either vanishing (cx to 2 ) or suppressed by 
e~ m l T , if the particle propagating in the loop has mass 



-m/T 

to. The explicit expression corresponding to diagram (a) 
reads 
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where the abbreviations u] n 

l/(e w "/ T — 1) have been used. The symbol A denotes 
the usual coupling constant of the linear a model, and 
Z-n- = i is the residue of the pion propagator at tree 
level §]. 

In the numerical calculations, the finite parts of the 
counterterms are determined such that f„ in vacuum re- 
mains unchanged by the loop corrections. In the chi- 
ral limit this has the consequence that Aero V ac van- 
ishes. Hence, we can replace oo in Eq. (^) by (ct)t=o- 
For further simplification we make use of the relations 



2A <Jq and / OT = <Jo /\/Z^, which follow from those 



given in Ref. || 
finally obtain 



in the limit = 0. In this way we 



I T 



(4) 



as required by chiral symmetry ||. 

Let us now turn to the temperature dependence of the 
p self-energy H p (qo,q), where qo and q denote the en- 
ergy and three-momentum of the p meson relative to the 
heat bath. In the following we will restrict ourselves to 
q = 0. In this case the self-energy tensor E^" has, as 
in vacuum, only two independent components, namely 
the four-dimensionally longitudinal and transverse self- 
energies, E p and In what follows we will concentrate 
on the latter and denote it by £ p for simplicity. Similar to 
Aoo, the self-energy E p can be decomposed into vacuum 
and medium contributions. The dominant medium con- 
tributions are generated by the diagrams (b) to (d). The 
two-pion loop, which is the dominant graph in vacuum, 
does not contribute to order T 2 , because the pions couple 
to the p meson in p-wave, resulting in additional powers 
of the thermal pion momentum. The explicit expressions 
for the three contributions lead to 



^plied ~ 2h 2 <7oAa med 

(c ) __ (2hi + 3h 2 )Z 7r 



J p med 



2tt 2 



d\k\k 2 



(5) 
(6) 



E 2L = ^/W-^=- 

/ U ai - 



,g 2 - {uo ai - uj-k) 2 q 2 - (w 0l + ^tt) 



(7) 



(To obtain real and imaginary part, qo must be replaced 
by qo+ie.) The contribution of diagram (d) takes a much 
simpler form in the chiral limit if, in addition, T -C \qo — 
m ai \. In this case, to order T 2 , the pion momentum k 
can be neglected in the second line of Eq. (0), i.e., = 
and u! ai = m ai , and the sum of the three contributions 
reads 



J p me 



<*(?<), o) = 



T 2 h lf j 2 
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(8) 



The change of the p mass is governed by the real part of 
the self-energy at qo w m p , but Eq. (@) implies in partic- 
ular £ pmec i(mp,0) = 0. Thus, to order T 2 , the p mass 
does not change at all, even if the condensate decreases, 
because the self-energy contribution oc Aaomed, diagram 
(b), is exactly canceled by the other diagrams, (c) and 
(d). 

For the a\ meson the situation is analogous. The domi- 
nant medium contributions to the a\ self-energy are given 
by the diagrams (e) to (g). In exactly the same way as 



described above for E p , one finds, to order T 2 : 



J (90,0) 



6/| ql - m 2 



(9) 



This shows that also the a% mass remains constant. The 
fact that the p and a\ masses do not change to order T 2 
was already found within the so-called gauged linear a 
model some time ago |Tc|| . However, within the gauged 
linear a model, this is less surprising, since the p mass 
is independent of oo at tree level, and diagram (b) does 
not exist. 

The results presented above are valid only in the chiral 
limit and at low temperatures. The more general case of 
TOjr 7^ and higher temperatures can be analyzed nu- 
merically. Then, of course, not only the diagrams shown 
in Fig. ^, but all one-loop self-energy diagrams have to 
be included. We will use two different parameter sets: 
one of those listed in Ref. || (parameter set A), and one 
with a very small pion mass (m» = 0.1 MeV) appropri- 
ate for the chiral limit. In both parameter sets we have 
chosen mo — 0, i.e., naively one would have expected 
" Brown- Rho scaling", as explained below Eq. (Q). 

As in Ref. || , we define the one-loop corrected p mass 
nip 1 ^ as the position of the maximum of the p spectral 
function, rather than as the zero of the real part of the 
p propagator. This implies that m p could in princi- 
ple change, even if E pme d vanishes at qo = m p . However, 
this subtlety does not change our conclusions for the low- 
temperature behavior of significantly, as can be seen 
in Fig. ^. In the chiral limit, to^ 1 ' stays almost constant 
at low T, whereas the condensate decreases rapidly ac- 
cording to Eq. ([|). At very high temperatures (T > 150 
MeV), of course, the low-temperature expansion is no 
longer valid, and also the p mass starts to change. In 
the case m w = 140 MeV, the main difference to the chi- 
ral limit is that all temperature effects are suppressed 
by e - '™"/ 7 , and therefore also the condensate remains 
almost constant below T w 50 MeV. At higher temper- 
atures, the cancellation of the different contributions to 
£pmed is not as perfect as in the chiral limit, resulting 

in a somewhat stronger temperature dependence of m p 
than in the case = 0. 

Our results corrobate earlier arguments (2| that in the 
chiral limit a change of the p or a\ mass of the order 
T 2 would violate the Eletsky-Ioffe mixing theorem and 
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FIG. 2: The temperature dependence of the condensate (a) 
and of the p mass m^ 1 ' in the chiral limit (solid lines) and 
for a realistic parameter set (dashed lines). All curves are 
normalized to their values at T = 0. 
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FIG. 3: Diagrams for the dominant medium contributions to 
E 7 in a strict loop expansion to one loop order. 



is generated by the diagrams shown in Fig. |^. (As in 
Ref. fa], we draw diagrams for the photon self-energy 
E 7 = — e 2 IIy rather than for Ily.) The contribution of 
the diagrams (h) to (j) can be expressed in terms of the 
p self-energy: 



n 



(Ih-j) 
V med 



(hi) 2 



(l 2 o 



n\n med ■ 



(15) 



is therefore forbidden by chiral symmetry. This mix- 
ing theorem, which relies only on the current algebra 
of chiral symmetry, states that to order T 2 the vec- 
tor and axial- vector current-current correlation functions 
Hv{qo,q) and U.a{QoiQ)i respectively, behave as follows 

mi-. 



Ily = (1 
n A = (1 - e)U A 



e)Hvvac + ^Avac , 
f ellv vac , 



(10) 

(11) 



withe = T 2 /(6/ 2 ). 

It should be noted that the mixing coefficient e is of 
order H, since the mixing is generated by thermal pion 
loops. Hence, if the correlators are calculated in a strict 
loop expansion, the mixing theorem cannot be fulfilled 



exactly. E.g., if Ily •* and 11^' ; denote the correlators to 
order h n (i.e., up to n loops), the medium contribution 



to n 



Y must be given by 



The diagrams (k) to (n), which can be evaluated with 
the same techniques as described for S p me d , give 
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and the sum of all contributions reads 
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(12) 



Since our Lagrangian is chirally symmetric, this relation 
is certainly fulfilled. At tree level the transverse parts of 
the correlators read 



<V) 



(fq 2 f 
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(13) 
(14) 



where / and g denote the 7/? and pirn coupling constants, 
respectively p|. The pion pole does not appear in Eq. 
( |l4| ) , because it contributes only to the longitudinal part 
of the axial-vector correlator. The last term in Eq. ( H ) 
can be understood as the contribution of the WWaa 
seagull vertex, if 11^ is computed via the M^-boson self- 



energy. Next we turn to the left-hand side of Eq. (12) 



i.e.. the medium contribution to IL;l . To order T 2 it 



which is exactly the result expected from Eq. (|l2|). 

However, in the model for HV and 11,4 we presented 
in Ref. ||, the correlators are not calculated in a strict 
one-loop approximation. Instead, the Dyson series is 
summed, allowing for a successful description of Ily in 
vacuum, but mixing all orders of Ti. Thus the dom- 
inant medium contributions to ITy are given by dia- 
grams similar to those shown in Fig. |^, but with all 
bare p propagators l/(q 2 — to 2 ) replaced by dressed ones, 
G pV ac(q 2 ) = l/[q 2 -m 2 p -T, pvac (q 2 )], and all point-like 7p 
vertices / replaced by loop-corrected ones, F lpvac (q 2 ). 

In this approximation one can see that, to order T 2 , 
the imaginary part of Ily (i.e., the spectral function) de- 
creases in the neighborhood of m p , but exactly at m p it 
remains constant, i.e., the strength of the p meson peak 
is reduced by reducing its width and not its height. In 



addition, similar to Eq. (20), a peak cx 8{q^ — m a ) is 



generated, but the coefficient is different from that in 
Eq. (|2^). For demonstration, we have computed ImllV 
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FIG. 4: The imaginary part of vector correlator, Im Ily (qo , 0) , 
in the chiral limit as function of energy, in vacuum (dashed 
line), as an extrapolation of the low-temperature behavior 
calculated for T = 25 MeV to T — 150 MeV (dash-dotted 
line), and at T = 150 MeV (solid line). 
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FIG. 5: The imaginary part of vector correlator, Im Ily (qo , 0) , 
for a realistic parameter set as function of energy, at various 
temperatures: T = (dashed line), T = 100 MeV (dash- 
dotted line), and T = 150 MeV (solid line). 



at T = 25 MeV and extrapolated the medium contribu- 
tion to T = 150 MeV as if there were no effects beyond 
order T 2 . The result can be seen in Fig. || together with 
Im Ily vac ■ The third curve in Fig. || shows Im Ily calcu- 
lated at T = 150 MeV. Obviously this looks completely 
different from the extrapolation of the low-temperature 
behavior. The peak at go = m a x is washed out by the 
thermal motion of the pions, which is so strong that the 
irai loop [diagram (d)] generates an imaginary part not 



only at g 



but also at g 



Another impor- 



tant contribution is given by the two-pion loop, although 
it is of the order T 4 as discussed before Eq. (|J). These 
two effects are the reasons for the broadening of the p- 
meson peak visible in Fig. ^. 

To conclude this letter, we display in Fig. || some re- 
sults for the vector correlator obtained with the realis- 
tic parameter set, i.e., = 140 MeV. In this case, of 
course, the mixing in the form of Eq. (|ll|) cannot occur, 



since all temperature effects are suppressed by e ~ m,T ^ T , 
and the axial correlator mixed into the vector correlator 
at energy go must be taken at energies below go — or 
above go + m w . Indeed, with increasing temperature, in 
addition to the broadening of the p-meson peak, a sharp 
peak at go = m ai — shows up in Imlly, which can be 
interpreted in this sense. Unlike the peak at go = m ai 
in the chiral limit, this peak is not smeared by the ther- 
mal motion of the pions, because it results from a kink 



in S p at the threshold go 



m,r, which persists to 



arbitrary temperatures. The existence of this threshold, 
and therefore of the peak in Imlly-, is of course an arti- 
fact of our approximation, in which the ai-mesons inside 
the loops have no width. In a more realistic model this 
reminiscence of the Eletsky-Ioffe mixing would reduce 
to a smooth enhancement of ImII\/ above the p-meson 
peak. In order to see this, an approximation scheme, 
in which also the particles inside the loops are dressed 
self-consistently, is desirable fl2]| . 

In order to draw quantitative conclusions for observ- 
ablcs which are directly accessible by heavy-ion experi- 
ments (e.g., dilepton rates), also the effects of finite three- 
momentum q, and, more importantly, of finite baryon 
density, must be included. The former is straightforward, 
but the inclusion of baryonic effects in a chirally symmet- 
ric way is a very difficult task and is postponed to future 
work. 
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